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Abstract 

An experimental measurement system was developed and implemented by the NASA Glenn Research Center in 
the 1990s to measure turbofan duct acoustic modes. The system is a continuously rotating radial microphone rake 
that is inserted into the duct. This Rotating Rake provides a complete map of the acoustic duct modes present in a 
ducted fan and has been used on a variety of test articles: from a low-speed, concept test rig, to a full-scale 
production turbofan engine. The Rotating Rake has been critical in developing and evaluating a number of noise 
reduction concepts as well as providing experimental databases for verification of several aero-acoustic codes. More 
detailed derivation of the unique Rotating Rake equations are presented in the appendix. 


Symbols 


a,b constants 

A duct area 

B number of fan blades 

C normalizing coefficient 
D duct diameter 

E duct profile function 

f frequency 

h fan fundamental harmonic 
J Bessel function - 1st kind 

k wave number 

M Mach number 

m circumferential mode order 

N number of pulses per fan rev 

n radial mode order 

p acoustic pressure 

P modal pressure amplitude 
P modal power amplitude 
Q weighting function 

R duct radius 

r radial co-ordinate 

t temporal co-ordinate 

s fan shaft harmonic 

V number of vanes 

x axial co-ordinate 

Y Bessel function - 2nd kind 

Greek 

A admittance 

r| reduced frequency 

4> mode propagation angle 
k eigenvalue 

p density 

ct duct hub-to-tip ratio 

9 circumferential co-ordinate 
a> fan shaft rotational speed 
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Q rake rotational speed 

u fan/rake speed ratio 

v acoustic particle velocity 

C, cut-off ratio 

Z Impedance 


Introduction 

The NASA Glenn and Langley Research Centers have been involved in several programs (Advanced Subsonic 
Technology - AST, Quiet Aircraft Technology - QAT) whose goals were the reduction in transport aircraft noise 
attributed to the engine. A component of engine noise is fan tone noise caused by rotor-stator, and other interaction, 
coupled to duct propagation. Interaction modes contribute to the farfield tone noise generated by turbofan engines 
(ref. 1). In order to reduce tone noise it is desirable to understand the generation of these interaction modes. Detailed 
knowledge of the modal structure may indicate engine modifications that can quiet the overall tone signature. It is 
difficult to determine the modal content solely from farfield directivity. Analytical modal description (refs. 2 and 3) 
of turbofan noise is well understood. Previous methods to experimentally measure duct modes used a concentric 
rake face or a reposition-able radial rake. These methods had drawbacks in that they either contaminated the 
acoustic field with their wakes or required lengthy, discontinuous data acquisition. 

In order to develop fan tone noise reduction techniques, a diagnostic tool was needed to measure the interaction 
modes of ducted fans. A concept for a continuously rotating microphone rake (refs. 4 and 5) was put forth that 
would be able to obtain all modal amplitudes in a reasonable data acquisition time. The key concept of the Rotating 
Rake is that by slowly rotating in a manner locked to the fan shaft, a Doppler shift is imparted to the duct spinning 
modes that is uniquely based on the mode physics. The resulting Rotating Rake system developed by the NASA 
Glenn Research Center has made several advances in the understanding of turbofan mode generation. In addition to 
directly evaluating a noise reduction technique, this diagnostic tool can provide an experimental database to which 
the results of aero-acoustic prediction or propagation codes can be benchmarked. This paper describes the theory 
and implementation of the Rotating Rake with historical examples. An appendix to the full paper, published as a 
NASA TM, provides derivation of equations unique to the Rotating Rake. 


Duct Mode Theory 

First, a summary of the duct mode theory with an application to mode propagation and generation is presented 
in order to understand the noise signature being measured. 

Cylindrical Wave Equation 

The propagation of sound in a circular duct with uniform flow is governed by the convective wave equation. In 
cylindrical co-ordinates this equation is written in non-dimensional form as: 

V 2 /?(r,©,x,Q = jj^ + M j p(r,&,x,t ) (1) 

with a general solution of 

P(0,r,x,t) = p mnf * E mn (k mn r)e^ +me±k ^ (2) 

where 

^ mn mn ^ mn t ^ m mn ^*) Qmn (^" mn Q] (3) 

is a solution to Bessel’s equation. The index m is defined as the circumferential order defining the number of 
pressure cycles in the circumferential direction. The index n defines the radial order, the number of pressure nodes 
in the radial direction. The nature of the radial rake indicates that of primary interest is the radial solution. 
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The constant, C mn , is a weighting factor that normalizes the arbitrary Bessel function profile to a desired 
physical property. (As an aside the selection of C nm will affect the calculated mode pressure amplitude but the mode 
power level will be a fixed value regardless of the choice.) The normalizing physical property chosen for Rotating 
Rake analysis due to a desire to compare results to a predictive tool is (ref. 6): 
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For selected physical cases (hardwall condition or no-flow) this result arises from the principle of orthogonality. 
For the general case, the E-tunctions based on the eigenvalue solution are not orthogonal. Flowever the above 
identity for C,„„ is also used for convenience in these cases. 

The eigenvalues are determined from the boundary conditions; in the standard case it is the hardwall condition 
that the acoustic pressure be normal to the duct inner and outer wall (a consequence of zero acoustic velocity at the 
wall). 
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( 5 ) 

( 6 ) 
(7) 


Note that while we are focusing on the radial solution, it is coupled to the circumferential and axial solutions 
through the mode order m and the eigenvalue K mn . 

The above boundary conditions show that real eigenvalues occur in the hardwall case and complex eigenvalues 
result from the softwall case; this is a concept that propagates throughout the analysis. The eigenvalues are 
determined from the boundary conditions. It is convenient to normalize the radial co-ordinate by the duct wall outer 
radius and define the hub-to-tip ratio as the inner wall radius normalized by the duct radius. Also, as convention the 
1st eigenvalue is called the zero radial, designated by n = 0. This allows for a short physical definition of a mode as 
the number of zero crossings of the pressure profile. The mode m = 0, n = 0, is uniform in pressure across the duct 
and is defined as the plane wave. Figure 1 shows a graphical representation of these modes with the regions of 
light/dark representing positive/negative acoustic pressure. Mathematically, an infinite number of radial modes 
(eigenvalue solutions) for each circumferential mode (also an infinite set) exist. 
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Figure 1. — Duct mode graphical representation. 


Figure 2. — Wave propagation in cylindrical duct. 


Propagation 

Although an infinite number of solutions exist, physical conditions allow only for a fixed subset of modes to 
propagate. The separated solution in the x-direction was given as: 
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( 8 ) 


The axial wave number is k mn . Examination of this equation gives three possible solutions: 

(') 2 7fRlc> pK mn \ k mn = a — > e ,kx = e lax => cyclic propgation in x - direction 

ii) 2 nfRtc = [3 k mn ; k mn = 0 (9) 

hi) 2 nfR / c < [3k mn \ k mn =ib — » e llcx =e'^ b ^ x =e~ bx => exponential decay in x - direction 

A real term when incorporated into the exponential term results in cyclical variation in the combined term. Thus 
a real wave number results in propagation in the x-direction. An imaginary wave number results in a pure real in the 
exponent, resulting in exponential decay. The special case of k mn = 0 defines the cut-off frequency. Physically, the 
3D waves have a helix propagation angle. Highly cut-on modes approximate a plane wave. Cut-on modes have an 
axial wavelength that is longer than would occur in free space at that frequency (shorter wave-number). Cut-off 
modes have ‘infinite’ axial wavelength (zero wave-number - decay). Figure 2 illustrates this propagation. The 
propagation angle (Ref. 7) is given by: 


cos(^) = 


—M ± 




1 + M. 1- 





for M = 0; 


where 


s=- 


2 jtfR 


cy/l- 


M ‘ 


= 2 7fRlK mn c for M = 0; (10) 


It is worth repeating the concept that a given (m,ri) mode will propagate only if the frequency is above the 
cut-off frequency. For frequencies only somewhat below cut-off, it can be shown that a 50 dB attenuation per duct 
diameter occurs. The eigenvalue, K„ m , increases with m or n. that is, the cut-off frequency increases with higher 
modes. This results in a sequential cut-on of modes from low order to high order as the frequency (or RPM) 
increases. For a given duct, the modes that propagate are from 0 to tn max ; the highest number of radial modes, from 
n = 0 to n = n max , (0,0) to ( 0 , n max ) will propagate at m= 0; each successive circumferential order will have fewer 
radial modes until at m max only the first radial mode will propagate (m max ,0). 
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Rotor-Stator Interaction 


The number of possible modes, while reduced from infinity, can still be substantial. A typical high-speed fan 
(shown later) may have up to 50 ///-orders with up to 12 radial modes resulting in 600 modes. Fortunately, the 
number of modes that need to be considered is a subset of these propagate-able modes. 

The classic paper by Tyler and Sofrin (ref. 8) presents the theory of fan-duct mode generation and propagation. This 
is seen by the kinematics between the rotor and stator. Figure 3 illustrates the progression of the pulse (which can be 
extrapolated to a continuous, cyclic profile). An illustrative rotor with B = 8, and V = 6 stators demonstrates that an 
interaction between a select rotor blade and stator vane progresses in the circumferential direction as the rotor turns. 
This pulse can be seen to travel in one complete revolution when the rotor has turned only 1/4 revolution. 

The rotor locked mode ( k = 0, m = B) spins at the shaft rotation speed. This mode, and higher order modes, can 
only propagate in a narrow annular duct if the blade tip speed corresponding to its spin rate is above Mach = 1.0. 
Lower order modes, which spin faster, may propagate if their spin rate results in a sonic blade tip speed. The critical 
tip Mach number is greater than 1.0 for non-narrow annular ducts. Sofrin demonstrated that generation of, and the 
spin rate of, circumferential spinning modes is governed by the following equations: 

generation: m = sB ± kV ; spin rate: O m = sBO. / m (11) 


Each circumferential mode, m, can have one or more radial modes, n. Mode propagation is dependent on the 
cut-off frequency, which is unique to each (m,n) mode. This frequency is dependent on geometric parameters and 
the eigenvalue of a combined Bessel function that is a solution to the cylindrical wave equation. Below the cut-off 
frequency the mode will decay exponentially. Above cut-off propagation occurs and acoustic power is transmitted 
down the duct and into the farfield. The Bessel function eigenvalue incorporates the duct geometry effects. The 
cut-off frequency, and cut-off ratio with Mach number effects is given by (ref. 9): 
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cut-off 


. CK mn 

nD 
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( 12 ) 


The cut-off ratio reduces to the familiar form when the eigenvalue is real (hardwall). The cut-off ratio rewritten 
in traditional turbomachinery form gives the ratio of the mode frequency to its cut-off frequency (for M = 0, 
hardwall case): 


ttsBOD _ sB (0/60) 

60/C nm C J cut-off 








Figure 3. — Interaction of rotating pressure pulse stationary object. 


(13) 
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Figure 4. — Duct mode generation and with propagation constraints. 


The primary method of generation of these spinning modes is the periodic interaction of the rotor wake 
impinging on the stator vanes as described. The stator vanes primary purpose is recovery of the rotational energy 
generated from the rotor (swirl). Other primary sources are inlet or exhaust duct struts, installed for structural 
support. These are typically the primary modes of interest in turbomachinery. Secondary methods of generation are 
caused by a change in rotor blade loading due to aerodynamic disturbances upstream of the rotor. If these 
disturbances are periodic with respect to the rotor, spinning modes may be generated. Discontinuities in the duct 
wall can generate unwanted modes. These modes, while generally undesirable, can contain significant acoustic 
levels. A summary of the preceding sections is given in figure 4. 


Rotating Rake Implementation 

Operating Principle 

The underlying principle of the Rotating Rake is the concept from the previous section that every 
circumferential mode spins at a different, discrete rate. The rake, rotating slowly, imparts a Doppler shift on the 
fan fundamental that is unique to each circumferential mode. This frequency separation is a linear function of the 
mode number and the rake rotation speed; the higher the mode number the greater the frequency shift with m = 0 
un-shifted (remaining at the harmonic frequency). This relationship, (derivation shown in the appendix) is: 

let — = — , free/ —>= oj(B - — ) (14) 

CO V V 

Figure 5 shows a spectrum from a single microphone with the rake stationary, compared to a spectrum for the 
same microphone with the rake rotating. The separation of the modes about the fundamental shaft order 16 is clearly 
seen. Cut-on background modes are 20 dB above the cut-off modes. The dominant rotor-stator interaction is 15 dB 
above these extraneous modes. 
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Figure 5. — Illustration of spectral Doppler shift due to rake rotation. 


Implementation 

The rake acquires unsteady pressure data as it spins using a series of radially distributed microphones or 
pressure transducers pointing at the fan. The rake is used in two locations — the inlet and the exhaust. The inlet rake 
rotates around the duct centerline inside the nacelle at the throat. At the nozzle, the rake rotates at the exit plane of 
the bypass duct. The rotational speed of the rake was required to be a specific whole number fraction, and in phase 
with the shaft. This means that when the fan has rotated u times, the Rotating Rake will have completed exactly one 
revolution within ±0.2°. The speed ratio, u, has been from 100 to 250. This speed and phase synchronization is 
necessary in order to resolve the Doppler induced frequency shift of the spinning acoustic modes. Consequently, the 
rake requires a very accurate drive and control system. This accuracy must also be repeated over the entire test 
period (5 to 10 min.) without build-up of phase error. The rake must follow the fan precisely even though the fan 
will wander in speed (5 to 10 rpm or so) throughout the test period. Consequently, the rake drive system must be 
slaved directly to the fan. 

The general method used for this precise following requirement is a programmable closed loop stepper or servo 
motor drive system. Both have been used in previous Rotating Rake applications. Commercial programmable motor 
control/drive systems are well suited for this type of application, easily meeting the design requirements if carefully 
selected. The rake motor drive system is an internal closed loop drive utilizing a programmed step/direction 
controller output driving a brushless servo or stepper motor; externally the motor drive system is an open loop 
system. The fan velocity is determined from a multi-per-rev physical sensor mounted on the shaft. A once-per-rev 
signal sets the fan position. Variations of this sensor have been an indexer mounted on the shaft, a gear tooth 
pick-up, or a blade -proximity sensor. The higher the number of signals per revolution, the better the rake system can 
track fan speed changes smoothly. This signal is fed into the manufacturer’s software for the motor control. A 
timing signal from the rake rotational body is used solely for position verification by monitoring on an oscilloscope. 
If the rake speed sensor resolution (pulses-per-rev) is equal to the speed ratio v, then when the rake is synchronized, 
the traces of the fan-one-per-rev and the rake u-per-rev will line up exactly, with no drift in time. Finally a rake 
body one-per-rev signal is used to provide information on the rake absolute position, which is important for 
determining absolute mode phase. 

The transducers’ time history signals must be brought across the rotating boundary. The method used is to bring 
the conditioned signal across by VHF/UHF telemetry. Again, carefully selected commercial units with custom 
design modifications have been found to be adequate to these requirements. The signal conditioning and telemetry 
transmitters must be mounted in the Rotating Rake drive housing. Power to supply these electronics can either be 
supplied by batteries in the rotating platform, or brought across by rails. The telemetry receivers must be diversity 
(two channels per frequency) type as the rotation will cause dropouts in the signal that must be picked up by a 
second channel. It has been found that close mounting of the antenna pairs is required. 

The output of the telemetry receivers is fed into a high-speed multi-channel data acquisition system. A primary 
requirement of the data acquisition system is the ability to synchronously sample the analog signal using the fan 
multi-per-rev signal, N. For analysis purposes, the ideal choice for N is 8B. This is because typically the 3rd 
harmonic ( 3B ) is the highest ‘frequency’ typically analyzed. To satisfy the Nyquist criteria, a sampling rate more 
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than twice that must be used. Further, since FFTs are more efficient when using a power-of-2, 8 B fits nicely. The 
signal must be synchronously acquired to the external fan shaft signal due to the time domain averaging employed to 
separate the tone signals from the broadband and flow noise. Think of the whole system as a strobe light. If the 
strobe light is set to a fixed (internal) rate, the rotating object may appear to drift. Connecting the strobe light to a 
timing signal causes the rotating object to ‘freeze’. Experience has shown that if any one of these synchronizing 
components is in error (drift, sample drop-out, noise) the data will be invalid. Best wiring practices are essential. 

Other requirements of the data system are the ability to record and store long continuous time histories 
(~5 min.) without breaking the trace into records. This requirement arises from the 25 continuous rake revolutions 
required to obtain good separation of the tones from broadband and flow noise. Ideally, each channel must have 
‘zero’ relative phase shift, or this must be accounted for in the post-processing. 

Finally, the number of transducers chosen must be greater than the highest number of radial modes that will be 
anticipated, usually the number of radial modes predicted to be cut-on at the m = 0 at the 3rd harmonic. A uniform 
radial distribution of the transducers along the rake has been chosen based upon experience. 

Rake Effects 

An obvious concern when taking inlet measurements would be the extraneous noise caused by the ingestion of 
the rake -wake contaminating the results. Indeed a stationary rake (V = 1) would generate every mode. However, in 
the reference frame of the rake, the rotation of the rake-wake and the relative speeds of the fan/rake cause all of the 
modes generated by the rake -wake rotor interaction to occur only at a single mode, that of the m = nB (fig. 5). This 
is also the so-called rotor mode that is cut-off below sonic tip speed. Therefore, the rake -wake interference is outside 
the typical modal analysis range unless the rotor tip speed is sonic. Conceptually, in the reference frame of the rake, 
it is as if the rake is stationary and the fan rotates at (f2-co) so that the rake-wake modes appear at that frequency. 

The mathematical verification of this intuitive observation is shown in the appendix. 

Both rakes have windscreens covering the transducers. The oscillating vortex shedding from the inlet rake body 
occurs at the measurement frequency of interest under certain conditions. A porous aluminum foam windscreen was 
attached along the trailing edge of the inlet rake covering the microphones to eliminate the vortex shedding noise 
from the rake body (the transducers are on the trailing edge of the rake) as seen in figure 6. In the exhaust duct, the 
rotor wake potential field convects downstream with a circumferential variation in a manner similar to acoustic 
modes. This can severely contaminate the microphone signals resulting in poor data. A tri-shield arrangement for the 
exhaust rake to lower interaction with the viscous wake was developed based on the concept that the convection is at 
the duct Mach number rather than the speed of sound resulting in a steeper propagation angle. The data quality was 
greatly improved with the addition of the shields as shown in figure 7 for a typical microphone. The convected wake 
levels at the non-propagating modes (m = -12, +16) is nearly eliminated, and the contamination of the propagating 
mode m = + 2 propagating mode is reduced. The resulting significant improvement in the radial profile is shown in 
figure 8. Finally, in the exhaust, the rake blockage must be taken into account. Typically for a high-pressure rise fan 
the nozzle area must be increased to account for this in order to make mode comparisons on a relevant basis. 


S= fd/V 

0.2 = (1 000 Hz X 0.02 ft)/(l 00 ft/sec) 


Si * 



FREQUENCY (Shaft Orders) 



Figure 6. — Reduction of strouhal-shedding 
due to installation of inlet rake shield. 


Figure 7. — Reduction of wake contamination 
due to installation of exhaust rake shields. 
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Figure 8. — Improvement in radial profile due to installation of exhaust rake shields. 


Data Reduction and Analysis 

The pressure time histories are acquired and recorded for each transducer as well as the timing signals through 
the data acquisition system described earlier. Experience has determined that 25 rake revolutions of data are 
required to obtain good signal-to-noise levels in the processed data. This time history is then time domain averaged 
(TDA) by taking record lengths 10 rake revolutions long and offsetting the next block by 1 rake revolution. This will 
typically result in 15 averages. Also this block length yields 10 bins in between modes when computing the FFT. 

The FFT of these TDA blocks are then reduced into complex magnitude for each Doppler shifted circumferential 
mode at each microphone radial location. A practical limit chosen for the modal decomposition is m = ±3 B, about 
each harmonic. The modal pressures are computed from the least-squares-fit of the radial Bessel functions to the 
actual pressure profile. A vector en'or is determined at each radial location from the difference between the actual 
and computed complex pressures divided by the maximum experimental pressure. The overall error of the solution 
is defined as the average of the radial mode vector errors. The discrete form of the pressure in the duct for a given 
frequency, m-order, at a specific radial location is: 


Pf 


,m,r = YjPf,m,n E mn (*,„„■ r) ! and in matrix form: [P r ] = [E w ][P n \ 


(15) 


n= 0 


applying least-squares analysis, inverting and solving for the modal amplitudes: 


(16) 


where 


N p 

N p 





X £ ™o (x ™ 0 r,) 

1=1 

1=1 

; K]= 

E m0( K m0 r l) 

Emo( K mO r Np ) 

(17) 

N p 

N p 


E,nN m . M XmiV max r l) 

■ • ■ E m N max ( K mN m3X r Np ). 


i= 1 

• • • ^ ^mA max ( K mmN max r i ) 

t=l 






P r is the vector measured pressures at a the microphone location and the solution vector, P„, is the modal 
coefficients; both complex. The size of the matrix is dependent on the number of microphones, and the highest 
radial mode being analyzed. Based on experience the numerical solution requires the curve fit to N max +1 radial 
modes for a hardwall solution. For a softwall solution the decay rate must be taken into consideration and may 
require a larger matrix. 
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Mode power is an important parameter in duct acoustics and provides a better description of the acoustic field, 
compared to mode pressure, especially when relating to the external acoustics. Power rather than pressure is the 
physical parameter that is conserved. Many wall modal techniques are unable to calculate modal power. The 
definition of power, and the final form of the computation used is: 


Power — » P = 



A 


(— + Uu)(p 0 u + puy, 
Po 


P = + nR (1 ] \P\* RE' 


Po c o 




1 + Mjl- 


I 


(18) 


This form mathematically illustrates the effect of cut-off ratio. For the hardwall case the cut-off ratio was shown 
to be a pure real number. If the cut-off ratio is greater than one, then the term in braces is also a real number, hence 
non-zero power. If the cut-off ratio is less than one, the numerator becomes a pure imaginary number due to the 
square root and power is identically zero. For the softwall, case the cut-off ratio is a complex number and non-zero 
power exists regardless of the cut-off ratio magnitude (power is the real component of the term in the braces). 

(z) C<h RE{bi\, P = 0 (0 \l\ < 1; REia+biV P*0 

Flardwall: real): „ Softwall: complex): „ (19) 

(ii) £> 1; RE{a}; P*0 (ii) \C\> l RE{a + bi}; P*0 

Other parameters calculated from the decomposed modal data that are presented are the pressure levels of the 
inner wall and outer wall. Cm, is a the normalizing coefficient discussed earlier. The function multiplier is the 
coefficient of the un-normalized Bessel function. The Sofrin coefficient is defined as the maximum pressure along a 
radial profile. The Power Level (PWL) is calculated, including cut-off ratio and Mach number effects, and 
summed to obtain the total mode power level. Table I shows typical output from the standard data reduction. All 
circumferential modes are analyzed in this manner and typically plotted in a 3-D ‘tombstone’ plot shown in figure 9. 
The base plane axes are and //-order, and the vertical value axis is the PWL in the (m,ri) mode. Along the wall of 
the /?7-order axis, the sum of all the radial modes provides the power in that circumferential mode. The sum of all the 
modes provides the PWL in the harmonic presented. A typical modal decomposition provides information as to the 
dominant modes present, usually those from the rotor stator interaction. Of secondary interest will be other modes, 
which may be due to inflow distortions or other geometric disturbances. Phase information, while tabulated is 
generally not plotted. 


TABLE I.— TYPICAL DATA OUTPUT 

Rate = MCFJnlet? 

Hub-to-Tip ratio = 0.000 

Corrected RR1 = 1800 Adjusted PPM = 1800 Duct Mach # = -0.115 
BPF Harmonic = 1 Circumferential Mode = 2 Padial Mode (max) - 1 


Location 

Measured 

Computed 

SfError 



Pa Deg 

Pa 

Deg 



0.943 

12.44 -124.2 

12.52 

-128.9 

8.23 1 


0.786 

11.50 -131.5 

11.30 

-132.2 

5.04 1 


0.628 

9.32 -132.2 

9.01 

-138.6 

8.58 1 


0.471 

5.93 -145.8 

6.16 

-146.1 

1.93 1 


0.314 

3.18 -156.2 

3.22 

-152.8 

1.55 1 


0.157 

1.26 172.5 

0.90 

-157.1 

5.37 1 


0.000 

0.28 45.7 

0.00 

0.0 

2.27 1 


Avg Vector ErrorSf = 3. 

14 






Mode Amplitudes 



Radial 

Cmn Real/Imag 


(Pa) Mag (dB) 

i Phase 


0 

2.720 -6.567 

-7. 082 

9.66 113. 

7 -132.8 


1 

3.342 -0 . 769 

0.458 

0.89 93. 

0 149.2 



Sofrin Coeff 

Inner Wall 

Outer Wall 

Radial 

Mag Phase Mag Phase tkg 

Phase 


0 

12.8 -132.8 

0.0 

0.0 12.8 

-132.8 


1 

1.5 149.2 

0.0 

0.0 0.9 

-30.8 


Mode Power 





Radial 

real/imag 

Cut-' 

off Pwr 

Decay 

PoTOr 

0 

1.780 0.000 

1.780 112.5 

0.0 

112.5 

1 

0.811 0.000 

0.! 

811 92.0 

17.1 

O.i 



Total Circum Mode Power = 

112.! 
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Achievements 


Several Rotating Rake systems have been custom built for three facilities. In order of complexity of the 
turbomachinery test article, these are (1) the Advanced Noise Control Fan, (2) various high speed scale model fan 
rigs in the NASA Glenn 9- by 15-foot wind tunnel, and (3) a full scale turbofan, the Honeywell TFE73 1-60. 

Advanced Noise Control Fan 

The Advanced Noise Control Fan (refs. 10 and 11) (ANCF) was designed in the early 90s and first checkouts 
occurred in 1996. This unique fan was designed as a test bed for evaluating noise reduction technologies, code 
verification/development, and to continually improve the Rotating Rake system. The ANCF is a 4-ft diameter ducted 
fan (fig. 10) with several unique features. The cantilever nature of the rotor support structure allows for the rotor 
alone to be run, eliminating the need for support struts. The 16-bladed rotor, in conjunction with variable count- and 
position stator vane sets, allows for a wide combination of modes to be generated. 

The ANCF has separate rotating rake systems in the inlet and exhaust ducts to allow for quick though not 
simultaneous measurements. The rotating rakes on the ANCF rotate at 1/1 00th of the nominal fan shaft speed of 
1800 rpm (i.e., 18 rpm) and contain 7 (inlet) or 6 (exhaust) microphones. A drive motor connected to a straight shaft 
and gear assembly is coupled to a ring gear. This ring gear is embedded into an Aluminum housing which forms part 
of the duct contour. These rake assemblies can be located between any of the discrete duct sections to allow 
measurements to be taken at several points along the duct axis. This has been extremely useful in validating 
propagation codes. Though the primary measurement system is the rotating rake, the ANCF has other aeroacoustic 
measurement systems (surface pressures, hot film, farfield) that allow for determination of the relationship of duct 
modes to the physics. 

The ANCF played a fundamental role in the NASA Advanced Subsonic Technology Active Noise Control 
(ANC) program. The evaluation of the modal effects of a tested ANC system provided by the rotating rake was 
crucial to understanding and improving the ANC algorithms. Several advances in the understanding of mode 
generation (ref. 12), propagation (ref. 13), and radiation (ref. 14) and verification of codes resulted from tests run on 
the ANCF. Other mode measurement systems (ref. 15) were verified on the ANCF providing confidence in those 
systems to be used in more demanding environments. The Hershel-Quinke (ref. 16) noise reduction concept was 
evaluated on the ANCF. These tests led to modifications of the H-Q devices for full-scale turbofan tests. A recent 
promising concept for fan noise reduction, wake filling through the use of trailing edge blowing (ref. 17) was 
evaluated on the ANCF and vetted for a more complex test on a high-speed fan. Advanced trailing edge blowing 
concepts are being currently investigated in the high-speed fan regime. 



Figure 10. — Rotating rakes on advanced noise control fan. 
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Twenty-Two-Inch Fan Drive Rig 

The 22 in. fan rig in the NASA Glenn 9- by 15-foot Anechoic Wind Tunnel is an air-turbine driven rig that can 
have a specifically designed fan installed and tested. Typically a test program will involve detailed performance and 
acoustic evaluations. The test section Mach number is generally held at 0.2. 

The first application of the rotating rake system (ref. 18) was in 1991 for the Advanced Ducted Propulser (ADP) 
fan model that was designed and built to test fan noise. This fan was an 18-bladed rotor 17.25 in. in diameter. The 
rake motor drive was connected to a large ring gear to which the rake body was mounted. The rake body/ring gear 
assembly rotated at l/250th of the fan speed (12,000 rpm) or 48 rpm. This assembly was fastened to the tunnel floor 
and ceiling via a support arm. The complete fan and rake assembly is shown in figure 1 1 . 

Configurations tested were a cut-on stator vane set (22 vanes) as well as a traditional cut-off set (40 vanes). 
Three inlet lengths were tested (L/Df = 0.21, 0.41, and 0.53). This test documented experimentally for the first time 
the effects of mode cut-on. A weak tendency for lower mode levels with shorter inlet length was noted. Also, it was 
discovered that very minor perturbations in the rotor casing treatment cause extraneous modes, resulting in higher 
acoustic farfield levels. This cause would have been difficult to identify from the farfield data. 

The next iteration of rotating rake mechanical design for the 22 in. fan rig was a conformal design to provide 
better mechanical stability and flow characteristics, and to allow for more flexible positioning of the measurement 
plane. The rake assembly was housed in a drum system that was mounted on parallel rails. The rake dram housing 
can be slid to various axial positions to accommodate different measurement planes; including all the way aft to 
enable exhaust duct mode measurements. The more robust structural mounting allowed for a faster rotational speed, 
l/200th of the fan. These features can be seen in figure 12. 

This system was used on the ADP in 1996. The purpose of the test (ref. 19) was to evaluate the potential for fan 
noise reduction by lowering fan tip speed. Acoustic modes were measured for two fans of equal pressure ratio. Fan 2 
had a tip speed 90 percent of Fan 1 depending on condition. Initial assumptions were that the lower tip speed fan 
would be quieter, but the mode measurements did not support this assumption. One complicating factor was that in 
order to maintain the performance Fan 2 required 5 1 stator vanes compared to Fan l’s 45. The different interaction 
modes generated (m = -9 vs. m = -15) at 2BPF may make comparisons difficult as the ANCF has shown that some 
interaction modes couple better to the duct. This test was the first to document the assumption that co-rotating 
modes (negative m orders) do not propagate well through the rotor. This transmission loss has a significant impact 
on inlet noise. 

The Active Noise Control program culminated in the test of a hybrid passive/active system (ref. 20) for tone 
reduction in a high-speed fan. The measurements made by the Rotating Rake were key to the design of the system. 
The baseline, passive only, active only, and hybrid active/passive systems were tested. The mode measurements 
showed significant reductions were achieved using the active control system, as well as a possible ‘synergistic’ 
benefit in that the reduction achieved by the combined system was greater than the sum of the reduction achieved by 
the systems individually. 



F igure 1 1 . — Rotating rake ring 
assembly on geared advanced ducted 
propulsor in the 9- by 15-foot wind 
tunnel. 



Figure 12. — Rotating rake dram 
assembly on 22 in. fan in the 
9- by 15-foot wind tunnel. 
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A Source Diagnostic Test (ref. 21) (SDT) was conducted in the 9- by 15-foot wind tunnel using the 22 in. fan 
drive rig by NASA Glenn as part of the AST. This fan was a 22-bladed fan that had sonic tip speed at 100 percent 
rpm. In one important phase of the test, the effect of stator sweep was investigated. The concept involved returning 
to a cut-on set of stators; sweep was introduced to mitigate the effect of cut-on BPF. The Rotating Rake 
measurements showed that compared to a standard cut-off set (26 vs. 54 Vanes), the swept vane set reduced the 
mode tone power significantly. This test also demonstrated the ability of the Rotating Rake to make valid mode 
measurements with the fan at sonic tip speed. Recall that at M tip > 1, the m = B mode cuts on. It was shown that in 
the inlet the rake wake contaminates the m = B mode. By analyzing mode measurements taken at fan speeds below 
the rotor mode cut-on, the magnitude of the rake wake contamination was estimated. When the mode measurements 
were analyzed above sonic tip speed, then increase in m = B mode PWL was greater than 10 dB. The physical 
assumption is made that the rake wake contamination would increase only slightly with rpm, but the rotor mode 
increases significantly. Therefore, the m = B mode measurements above sonic fan tip speed in the inlet are assumed 
to be valid. 

The rotating rake was also used to evaluate the acoustic benefit to a forward swept fan (ref. 22), the Quiet High 
Speed Fan (QHSF), compared to a baseline fan representing the Honeywell TFE731-60 engine. Being a 
representative fan model, support stmts were included in the exhaust duct, in addition to the stator vanes. This test 
also presented an extension of the Rotating Rake reduction methodology to the analysis Multiple Pure Tones 
(MPTs) or “buzz-saw” noise. MPTs occur at shaft orders below BPF with the interesting characteristic that the 
m-order generated at a given shaft order is numerically, the shaft order. The mode measurements showed that the 
QHSF was significantly quieter, but that the reduction was primarily from a lowering of the rotor -strut interaction 
tones rather than from a reduction in the MPTs that had been originally been the design goal. The Rotating Rake 
measurements suggested that the reduction mechanism may have been that the stator was shielding the stmt from the 
wake of the rotor; though this is unconfirmed. 

Follow up tests to the SDT, comparing different rotor sets; and the QHSF, measuring the effects of vane 
clocking, have been done, but not yet reported as of the date of this publication. 

TFE-731-60 

The most challenging and perhaps most significant application (ref. 23) was on a Honeywell full-scale turbofan 
engine, the TFE73 1-60. The TFE73 1-60 engine is currently certified on the Dassault Falcon 900EX. The TFE73 1- 
20/40/60 engine family also powers the Lear 45, IAI Astra SPX, and the Dassault Falcon 50EX. This test was 
accomplished under the Advanced Subsonic Technology program sponsorship of the Engine Validation of Noise 
Reduction Concepts (EVNRC). The basic engine parameters are listed in table II. The schematic of the fan flow- 
path in figure 13 shows multiple sources of interaction modes. The engine was tested at the Honeywell engine test 
stand in San Tan, Arizona. Photographs of the inlet and exhaust duct rake assemblies with the engine on the test 
stand are shown in figure 14. 


TABLE II.— TFE73 1-60 ENGINE PARAMETERS 


Basic Engine Weight 

988 lb 

Takeoff, Sea Level, Static Thrust 

5000 lb 

Airflow 

187 lb/sec 

Cruise Bypass Ratio 

3.9 

Cycle Pressure Ratio 

22 

Geared Fan Pressure Ratio 

1.70 
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Figure 14. — Rotating rake assembly on Floneywell TFE731-60. 


The rotor-stmt interaction modes at BPF were identified as the dominant modes but were a minor contribution 
to overall PWL, due to the strength of the extraneous (background) modes. Rotor transmission losses were shown to 
result in lower co-rotating modes in the inlet compared to the exhaust. The 2BPF rotor-vane interaction mode was 
clearly identified in the inlet and exhaust, and observed to dominate. Repeatability of the measurements was shown 
to be excellent by comparing results from two separate test builds. An attempt to measure the rotor-locked noise was 
inconclusive due to the weakness of the rotor mode. The MPT sub-harmonics (“buzz-saw”) were seen to be 
dominated by the shaft-orders just under ABPF. 


Future Extensions 

An evaluation of the Trailing Edge Blowing concept on the 22 in. Fan Drive Rig is underway in the NASA 
Glenn 9- by 15-Foot Wind Tunnel. Currently, a Rotating Rake system is being designed for a Floneywell Tech7000 
Turbofan (for the NASA sponsored Engine Validation Noise Reduction Technologies (EVNERT) program. The 
rotating rake will be used to measure the full-scale version of the QF1SF and to provide detailed design information 
for FI-Q tubes. The analysis technique is currently being extended to measure and reduce duct modes over treated 
wall sections. The method has already been verified for treated inlet ducts (a = 0), and current testing for annular 
treated sections is underway. These results will be reported in a future paper. A significant need is to develop 
analysis techniques to properly measure duct modes in sheared flows. 
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Appendix 


This appendix outlines the derivations of the equations presented in the main body of the paper. General 
assumptions are uniform Mach number in a constant area cylindrical or annular duct. 


a,b constants 

A duct area 

B number of fan blades 

C normalizing coefficient 

D duct diameter 

E duct profile function 

f frequency 

G mode power Mach effect 

h fan fundamental harmonic 

J Bessel function - 1st kind 

k wavenumber 

M Mach number 

m circumferential mode order 


Symbols 

N number of pulses per fan rev 

n radial mode order 

p acoustic pressure 

P modal pressure amplitude 

P modal power amplitude 

Q weighting function 

R duct radius 

r radial co-ordinate 

t temporal co-ordinate 

s fan shaft harmonic 

V number of vanes 

x axial co-ordinate 

Y Bessel function - 2nd kind 


Greek 

A admittance 

s acoustic particle displacement 

r) reduced frequency 

<t> mode propagation angle 

k eigenvalue 

p density 

a duct hub-to-tip ratio 

0 circumferential co-ordinate 

co fan shaft rotational speed 

Q rake rotational speed 

T Bessel function combination 

A duct mode amplitude 

u fan/rake speed ratio 

v acoustic particle velocity 

C, cut-off ratio 

Z Impedance 


A. Solution to Wave Equation with Boundary Conditions Representative in Turbofans 


The non-dimensional form of the convective wave equation in cylindrical co-ordinates with axially uniform mean 
flow in a constant area duct is written (ref. 1) as: 

V 2 p(r,0,x,O = jj^ + M -^-j p(r,@,x,t ); 

* * * * * * 

• ,7 , , u x r x t p , / R 

with the non - dimensionalizations M = — r ; x = — ; r = — ; r = — ; t = ; p = — t = — q— 

c R R R V. pc 2 c 

n / C • n n 


A. 1 


This form of the equation is known for it’s separable solutions: 

p{6, r,x,t)= p* E{Kr)e ia * t+me±kx) 
A cyclic temporal variation is superimposed upon the entire solution: 


P t (t)=e llnfi 


A.2 

A.3 


The form of the circumferential solution is deduced from the observation that the pressure must be single-valued 
azimuthally: 

P s (0) = P g (d +2 k)\ results in P g (<9) = cos (m 6) A. 4 

The variable m is defined as the circumferential mode order and indicates the number of pressure cycles in the 
circumferential direction. 
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The solution in the axial direction, where k is defined as the axial wave number, is given as: 


p A x )= 


—ik r ; 


K = L 

>1 (l-M')L 


-M± Jl-(1-M 2 )(-) 2 


rj= fRIc 


A. 5 


Since the text clearly illustrated that the circumferential solution arises from the rotation of the rake inducing a 
Doppler-shift that separates the circumferential modes. The nature of the radial rake indicates that the primary 
interest is the radial solution. The form of the radial solution is the well-known Bessel’s Equation: 


d 2 P(r) 1 dP(r) 

— + — + 

dr" r dr 


2 


n 


1-M — 

V rj) 



\ 1 ) 


( 

m 

~ 

v 


■\ 


J 


lP(r) = 0 


A. 6 


The solution to this ordinary differential Equation is an eigenvalue problem with an infinite number of solutions: 

P{r) = aJ(K n r)+bY{K n r)\ a, b arbitrary, n = 0,1,2,... A.7 

where J and Y are the Bessel’s functions of the 1st and 2nd kind, of order n. The index n is the radial mode order, 
which physically indicates the number of pressure nodes in the radial profile. The pressure profile is a superposition 
of cirumferential and radial modes. Each term in this superposition contains a radial duct profile function written as: 

E mn <r) = c mn x ¥ mn O) = C mn [J m (K mn r) + Q mn Y m (K mn r)J n = 0,1,2,... A. 8 


with subscripts m for the circumferential mode order, and n for the radial mode order. 

Since the Bessel functions can vary substantially in form, by convention a=C, b=C*Q is chosen where C is a 
weighting factor related to the integration of 'A and Q is termed as a second eigenvalue, though technically it is not. 
K"is an eignevalue found without reference to Q as shown in Eq. B.7 or Eq. B.8 solves the proper eigenfunction P(r). 
Q is from the boundary conditions and determines the proper radial shape of the eignefunction. 


The general boundary conditions are determined from the relationship between the non-dimensional particle velocity 
and pressure at the duct walls: 


f k V 

1-M — 

V J 


£ = 


dp 

dr 


p Z 1 

( 7% />o c o ^ 


A. 9 


Where £ is the particle displacment in the flow positive in the radial direction, and is the radial particle displacment 
at the wall, positive into the wall. When the boundary conditions are applied, this becomes: 


dP{ 1) 
dr 


= -iV A ow 


1-M — 

v n) 


P( i); 


dPia) 

dr 


ir/A h 


' k 

1-M ^ 

>1J 


P(c J) 


A. 10 


The boundary condition for the hardwall case is that the pressure be normal to the duct inner and outer wall (a 
consequence of the zero acoustic velocity at the wall). Or, for a hardwall duct the inner and outer wall admittances 
are 0 so Eq. A. 10 reduces to: 


dP(\) _ dP(a) _ Q 
dr dr 
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B. Eigenvalue Solutions 


The solution of the eigenvalue problem is based on the general boundary conditions from the continuity of particle 
displacement at the wall (Eq. A. 10) that are repeated here: 


dPj 1) 
dr 


~ ir lK w 


1 — M — 

v n - 


pay, 


dP(a) 

dr 


i n A lw 


1-M — 

v n , 


P{C7) 


B.l 


where. P mu {r ) E m ^P mn (^r) C mn \J m Qc mn r) + Q mn Y m (j^ mn r)] 


B.2 


Taking the derivative of this equation with respect to r, and applying the recursive definition of the Bessel Function 
derivative ( mJ m (ic mn ) - K mn J’ m (/c rnn ) = K mn J m+l (K mn ) ) results in: 


dE mn (r) 


dr 


d Cmn \_J m ( K mid ) + Qnm Y m ( y m;i r )] _ ^ r f (tr r\ + O Y’ (k r'll 

, mn^ mn L * 7 m v*- mn ' ) ' *£mn 1 m vA- mn ' ) J 


Y \ 



y \ 

n 

m 

J m fa mn ^*) — ^ m+ 1 fa mn ^*) 

+ Qmn 

m 

famn — ^m+1 famn ^*) j 

_U mn r ) 

U mn r ) 


B.3 


Substituting Eq. B.3 into the left side of Eq. B.l, and Eq. B.2 into the right side, canceling out the coefficient C m 
results at the outer wall: 


f ^ 

m 


r m ) d, n + 1 ) 


V /V mn J 


= -iijA 0 


+ Qmn 

y \ 

m 

pm fa mn ) ^m +1 fa mn ) 


A K mn ) 



' k 

1-M — 

\ V J 


\d m (Kmn ) Qmn Y m )] 


and at the inner wall: 


ff ( > 



Y N 

I 

tn 

J rn ( K mn a ) ~ J m + 1 

+ 2m« 

YYl 

Y m( K mn Cr )- Y m + Mmn G ) j 


\ K mn a ) 


= +iV A m 


f k V 

1-M — 


[J m (K mn v) + Q mn ^/nn^] 


v n J 

Solving both Equations for Q mn and setting the result equal to one another yields: 



J m famn ) ^ m + 1 ) 

1 

f-> 


m + M 0W 1-M — 

^m+1 ) 

v nJ 



m - 1 

iT vi 

a ' 

1-M — 

v h J 

2 

^ m famn ^ m+\ famn 

m 

( k v" 

2" 


— 

a 

1-M — 

v q / 


famn ^m+1 famn 


B.4 


B.5 


B.6 


This equation can be solved iteratively by a number of means to obtain k , then substitute into Eq. B.4 or Eq. B.5 to 
obtain Q. For the hardwall boundary conditions the admitances, A iw and A ow , in Eq. B.6 are set to zero: 


m (P- mn 2 ^ mn ^ m +1 ^ mn 1 

m Y m ( K mn ) - K mn Y m+1 (K 

mn ) 


d m ( K mn G ) ~~ K <mJ 1 ( K mn a ) 


'Pm f*'™ 17 ) - ^mnfa m + 1 (K mn <j) 


B.7 
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The profile for a cylindrical duct (ct= 0) requires that the radial solution be finite at the centerline. To Q=0 to Thus 
the radial duct function, E, reduces to: 


^/nn C mn E mn (t) m^mn^ ^ 

The cylindrical duct eigenvalue solutions for softwall boundary conditions are found by noting A iw in Eq. B. 1 and 
Q mn in Eq. B.4 are zero. This results in: 


L „mv (1); am __ 0 

dr V 1 / dr 


J m + 1 ( Emn ) 


J m (?mn ) 


-m- irjA ow \ 1 - Af — =0 


Finally, the eigenvalue equation for the cylidrical hardwall duct is found by setting A ow in Eq. B. 10 to zero: 

„ J m+\( K mn ) _ m _ A 

K mn , . , m- U 

m v*- m/i ) 


C. Derivation of the Weighting Factor, C m „ 

Since the profile of the Bessel functions can vary significantly with order, it is desirable to normalize the duct profile 
by a physical constant in order that meaningful comparisons of the mode amplitude can be made. The weighting 
factor can take several forms based on the desired physical reasoning. One form chosen in early modal analysis 
(ref. 2) was acoustic intensity based normalization by setting the integration across the duct to unity. 

intensity: I p 2 dA = \ C.l 

A 

An alternative form is to normalize by relating to power. In order to match the methodology of standard CAA codes 
(ref. 3) the normalization chosen is power based so that: 

power: J p~dA = A C.2 

A 

From the section A, the radial profile is given by : 

(r) = J m i K mn Y ) + Qmn Y ,„ (*mn * = 0,1,2,... C.3 


integrating, 


\ 2 J J 1 C 2 m A (K mn r)rdrd® = 2v(l - a 2 ) 

" 0 * (7 

C mn J X * (' K mn r) >'dr = (l - CJ 2 ) 


a-<r z ) 


0 Kmn r )rdr 


Under certain conditions the functions are orthogonal. The orthogonality property of the Bessel Function (ref. 4) is: 


if j * l 




If ttfMdt, ifj = l ; 


provided the following two conditions hold : 

(1) X[ is a real zero of /? 1 X^ V+1 (Kb) - h 2 (Kb) = 0 

(2) there exisits k\ & k 2 ^ 0 such that (K/a)- k 2 t, v (K/a) = 0 
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Applying the hardwall boundary conditions from Eq. A. 1 1 and the recursive definition of the Bessel function shows 
that (1) & (2) are satisfied with the eigenvalues k. Re-arranging the impedance boundary conditions Eq. A. 10: 


d(E mn (K mn r)) = _. i(1 _ M {k ^ r) ^ {k ^ r) + iA{l _ M {Kmn r) = 0 


^'m 2 , 

>1 


dr // // 

At the outer wall (b=l) this equation takes on the form of condition (1): 


C.6 


f K mn^m+l,n ( K mn ) m ^mn ( K »w )1 ^ ) ^mn^mn ) ^ 

L J r| 


K mn^m+l,n ( K mn ) 


x \2 


m + iA ow (\-M —) 


^mn ( k «hi ) ^ 


C.7 


Likewise, at the inner wall (a=o) the boundary conditions take on the form of condition (2): 


~ K mn E m+\,n ( K mn CT ) + E mn C K mn <J ) 

a 


- iA iw (1 - M —) 2 E mn ( K mn a ) = 0 

n 


K mn E m+l,n ( K "/n/i Cr ) 


— -iA iw (l-M^) 2 
a ?i J 


£ m »<Am„c)= 0 


C.8 


It is recognized that this is not strictly correct because for the general case the K mn is not a real eigenvalue, but may 
be complex. Nonetheless, the form of the solution will be used in order to be consistent in the normalization. Future 
work may investigate the derivation of a rigorously correct solution of the normalizing constant. 

The general solution to the integral is then: 


1 


( 2 ) 



| ¥« 2 ( K mn r ) rAr ~ 2 r ~ 


m 

2 2 

^ mn mn ( K mn r ) 


a 

. 

V *^mn ' ^ 

_ 



d-cr) 




f 2 



m 

2 

^ mn ( K mn ) ^ mn ( K mn ) 


V K mn ) 





K mn ® ) 


C.10 


Applying the general annular duct, non-zero impedance boundary conditions: 


dP( 1) 
dr 


= ~i>lA„ 




n ) 


dr 


i n A h 


' k ' 

d J 


\P((T)- 


cEV n M,nn r ) 

dr 


= —i rjA 


\ — (K mn r) 

n . 


K m ,Xn (K mn ) = -ivA 0 


' k ' 

1 -M^L 

U 


Wm ( K mn )i K ,„Xm ( K mn°) = ir lX 


\-M —\V m {K mn cr) 

n. 


u v ( K \ — — 

A m v /v mn / 


ir)A 0 


1 -M^L 

\ n j 


P^m mn ) ’ ^m (^* mn 


ir/A h 


l-M — 
>1 J 


WmtemnO) 


C.ll 
C.12 
C.13 
C. 14 
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Substituting the re-arranged boundary conditions equations into the definition of C m „ derivative: 


Cln 2(1 -a 2 ) 


.2 


1- 


m 


u/2 / 

x mn V inn ) 


f -n A f 
n A ow 


f 

2 

\ 

a 2 

m 


7~ 



_v 

^mn 

/ 


^ mn ( 


V 


M, 


1_M-L ^- m (K mn ) 
11 


1-M — | 'F m 2 (K m „a) 

il 


For the hardwall condition A=0; the definition of C mn reduces to: 
1 1 


( 2 ^ 

1 111 

u/2 A- )- 

1 mn \ /v mn ) 

f 

a 2 - 

2 ^ 

m 

' V ln( K m„v) 

2 

2 

_V K-mn J 


V 

^ mn J 



cL 2(1- a 2 ) 


For a cylindrical duct, equation C.15 simplifies to (impedance condition): 

with (7=0; X,u, ( K mn ) = J,n i K mn )> ^mn (0) = hi (0) = 0; III * 0 

\ 2 s , n2 


C, 


1 _ 1 
2~ _ ~2 


1- 


m 


7 2 (k 1- 

J mn 1 mn J 


Mo 


\ ^ mn j 


1-M — 
il . 


^ mhmnl 


and for the hardwall condition, this further simplifies to: 


1 1 n m ” 2 , . 

_? ^ ? 'J m ^ mn ) 

Ct„ 2 Km „ 


D. Doppler-Shift of Duct Modes in Rotating Reference Frame 


oo m=+M 

The pressure at a fixed radial and axial location in the duct from all modes is: p(®,t ) = I 

5=1 m=-M 


co m=+M 


If the location is slowly rotating in the circumferential direction: 0 = Qr; p(@,t ) = Z 2> 


s g-i(sco-mQ.)t 


5=1 m=—M 


For a typical application the rake rotates at a fraction of the fan speed: = — ®; s= B 

v. 

Substituting Eq. D.3 into Eq. D.2 it is seen that the cyclical part becomes: 

(jco-raO) = = 0) ^-( 7 ^/)) 

The mode frequency is then: 


/« =»(*-('%)) 
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C.15 

C.16 

C.17 

C. 18 

C.19 

D.l 

D. 2 
D.3 

D.4 

D.5 



Thus the frequency is shifted from the fundamental by an amount proportional to the mode number, and inversely 
proportional to the speed ratio. Note that for a rake that rotates in the same direction as the rotor, co-rotating modes 
(+/;;) will be reduced in frequency and counter-rotating modes (-/;/) will be shifted higher in frequency. The plane 
wave ( m=0 ) is at the fan harmonic. 

E. Non-Interference of Rake Wake 

The wake of the rake is ingested by the rotor and contaminates the acoustic field. Fortunately, due to the rotating 
frame of reference all of the contamination occurs at a single Doppler-shifted frequency. A physical explanation was 


given in the text; here is the analytical proof: 

The pressure due to a fixed distortion: p(®,t ) = A s m e ls±l)0) E.l 

If this distortion is slowly rotating: p(®,t ) = E.2 

The cyclical portion is: ( sco±lQ)t - (s± 1)9 = ( sco±lQ)t - ( s±l)Qt = scot ± IQt- sQt + IQt = s(co-Q)t E.3 

so the frequency (with s=B ) is: f rake = s(co - Q) = cos (l -Q/co) = coB( 1 - y') E.4 

from the previous section note that for m=B : f m= g = co j B - (^/) j = wB (l - (J/) j E.5 


In an unpublished letter, dated 20 November 1981, (this was to be an appendix to an earlier report (ref. 5), Sofrin 
demonstrated the non-interference this way (original text and notation): 

“Consider a single mode generated by the rotor cutting the wake of a radial probe positioned at 9=0: 

Pnm (0; 6, t) = Re A mn exp i (m 9 - nBQ.t) ( 1 ) 

If the wake is moved to 9=p, then the pressure at (9+[:S) and time (t+ p/Q) will be the same as given by (1), or 

p mn (A6+p,t+%)=p mn m,t) (2) 

Consequently, the pressure at 9 and t due to the wake at 9-p is 

Prnn (P; 0- 0 = Pnm (0; 9 - P, t - %) = Re A mn exp i m( 9 - P) - nBQ(t - (3) 

The complete pressure, due to all modes and frequencies is then 

Pnm (A 0, 0 = X d - ? ) = Re X A mn eX P ' [ m ~ P) ~ nBQ - f + nB P)\ ( 4 ) 

m,n m,n 

But the microphone senses the pressure only at the position 9=p, (or 9=P+const) so the microphone signal p mic is: 

Pmic = Pmn t) = Re^ A mn exp i \nBfS - nBOt] (5) 

m,n 
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Since the exponential does not depend on m , this can be written 






Pmic iVV / / 1 * 1 mn 

=1 V —oo J 


exp i (nBj3 - nBflt ) 


00 

or, with YjA„ m = A n , 

-00 

00 

Pmic = RelA, exp; iiB(fj-Dj) 

n = 1 


( 6 ) 


From Eq. (6) it is seen that the microphone signal, at any frequency n/i/2 behaves exactly as if it were sensing the 
direct rotor field which has nB lobes, spinning at £2. This signal will undergo a 2k phase shift when the probe angle 
(3 changes by 27i/p. 

It should now become clear on reflection that no matter how complicated the wake-generated acoustic field, that 
after the probe has turned through 1 blade gap, the pressure sensed by the microphone is indistinguishable in 
amplitude and phase from its initial state, thus providing intuitive support to the analytical results.” 


F. Least Squares Fit 

The modal solution is a least-squares curve fit to the measured pressure using the analytical radial pressure profile 
solutions as the basis functions. If the pressure has been separated by frequency and circumferential mode order, the 
pressure for at given radial and axial location is then: 

00 

Pmrf = YjPmnf E mnf ( K mnf r ) F - 1 

n= 0 

N max 

Limiting to the maximum number of radial modes: p mr f = /, p m „ f E mnf (K mnJ -r) F.2 

n= 0 


N max 

The difference between the measured and calculated pressure: 8 = p m „f E „mf ( K mn f r ) ~ Pmrf 

n= 0 


F.3 


A weighting function is introduced to allow for weighting of errors based on area or other physical means. (In 
practice the weighting function W t has the effect of bad coding a microphone by setting W t = 0 for a suspicious 
channel, 1 for all others.) 

N ma x 

^ 1 Pmnf^mnf ^ mnf Pmrjf 
- n = 0 


8i=W, 


F.4 


The square of the error is: 


sf = W, 


iVmax 


^ . Pmrf E mnf mnf ) P 


"’’'if 


F.5 


JV„ N„ 


Summing over the number of points p\ S 2 = w\ 


i= 1 i=l l L n=0 


N max 


P mnf E mnf C K mnf fi) Pmrp 


F.6 
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The concept of least-square method is to minimize the error squared, by taking the differential with respect to each 
P m „/ and setting the result to zero. Which taking the series terms out of the summation and extending to the radial 
modes: 


d5 2 

^Pmn\f 

db 2 

^Pmri2f 


Np 

2 Z^ 2 ( E mn x f^mn x f r i)) 
2=1 

N max 

(^Pmnf^mnf ) — Pmr\f ) 

22=0 

Np . 

A^max 

2£ W r ( E mn 2 f( K mn 2 f r i)j 

I i^Pmnf^mnf (^ 22222 / r i ) — Pmrif ) 


i = 1 


72=0 


F.7 


Setting the derivative to zero and splitting the summation: 


Np N max Np 

^ — ( E mvf fernyf f i ) ) Pmnf E mvf i^mnf r i ) ~ Z , [ E mvf ( K mvf r i )) [Pmrjf ) 

i=l n= 0 7=1 

N max Np Np 

Pmnf E mvf (. K mvf r i ^ E mnf 0^-rnnf r i ) — ( E mxf ^ myf r i )) ( Pm/}/ ) 

77 =0 7=1 7=1 


Which can be written in matrix form for a fixed frequency: 

[F m ][PnME vr ][P r ] 


[ / »]=[ / ;»r'K][/ j f ] 


F.9 


where: [F vn ] = 


n p 

Z 

7=1 


^1 E m o(K m O r i) 


N„ 


Z, W l E mN max v r i ^ E m0( K m0 r i ) 


2=1 


N P 

Z i E m0 ^ E mN maK ( K 777A r max r i ) 

7=1 P 


Np 2 2 

Z ^A7„ E , nNrr.™ ^ KmmN max r ‘ 1 


7=1 


y p mi v max 


F.10 


K] = 


^1 E m(>( K mO r l) 


W Np E mo( K mO r Np) 


l W l ^ miV max ^" lW max ^ - ^ Np E mN max (*" n)7V max r Wp )J 


measured pressures: [P,.] modal amplitudes (solution): [P„ ] 

note that the matrix Equation dimensions are: [N max ] ~ [N max , N raax ][N max ,N meas ][N raoas xl] 
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G. Derivation of Power 


The basic definition of power is the integration of the acoustic flux across a radial plane of the duct: 

P= J (l x )dA; 

A 

For uniform axial flow it is straight-forward to show (* represents the complex conjugate.): 

I x = (— + Uu)(p Q u + pU)\ p = c\p, M =U /c 0 — > I x = (l + M 2 )pu +——pp + p 0 c 0 Muu 
Po PoF o 

Inserting Eq. G.2 into Eq. G.l results in: 


P 

x mns 


2n 1 n2 /i _ 2 \ OO 00 00 

f J (l ax )rdrdQ = - ^ > Z I I G mns \P fj 


0 a 


p 0 u 


2 _ 

mns I 9 


m =— oo n = 1 s=- oo 


where G„ 


and A„ 


M~ -(1 + M )A mns + A mns 


sBQ. 

U 


-+y« 


• v = _r 2 k 2 

> I mns , I ' ' I mn 

G) 


G.l 


G.2 


G.3 


The derivation of this is shown in detail in reference 3. It is recognized that some error is introduced since the above 
result relies on orthogonality to eliminate cross terms. Summation of the +/- terms yields: 

G mns = M " - (1 + M - ). A mns + A mns A mnsmns 


2 2 * * 

* 

A- mn(-s ) A mns 

G mn(-s ) = M ” - (1 + M 2 )(A „ ms ) + A mm (A mns ) 


G.4 


G mns + G -mn(-s) - 2M ? - (1 + M “)( A mns + A,„„ s ) + 2A mns A mm 


G.5 


Substituting the expression for A mns into Eq. G.5, after manipulation: 


.. + M 2 j3 4 (sBn/U)k mns . 
( sBn/c 0 ±Mk mns f 


G.6 


Traditionally, this has been reduced assuming either pure real or pure imaginary eigenvalues to the form: 

G mns = . f or a rea | e ig enva | ue G mns = 0; for an imaginary eigenvalue. G.5 

(sBn/c 0 ±Mk mns ) 

A slightly different form will arise if the assumption that complex eigenvalues can occur: 
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M 2 /3 4 Re-] ~ (+k mns ) 


fl ___ 

^ mns 


CO 


G.6 


a> 


±Mk„ 


^2 j 2 

The final version for the rotating rake data analysis is: P = + \ 0 \ RE{ 

A) c o ' 


1- 


C 


G.7 


l±M D &-f 


F. Comment 

Derivation of the equations unique to the Rotating Rake system have been presented as used for several 
applications. The data reduction based on these equations has been validated by experience and comparison to other 
methodologies. In order to keep the solutions manageable a significant simplification assumption was made, that the 
orthogonality solutions for various integrals applies. For certain boundary conditions (namely flow in treated ducts) 
this is not rigorously correct, as it assumes product terms involving non-identical modes due not contribute. 
Flowever, this assumption provides a solution that can be meaningfully analyzed. Future extensions of this analysis 
should include the solution complete with cross terms that were not included herein. 
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